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ABSTRACT 

Thin accretion discs around massive compact objects can support slow pressure modes 
of oscillations in the linear regime that have azimuthal wavenumber m = 1. We con- 
sider finite, flat discs composed of barotropic fluid for various surface density profiles 
and demonstrate-through WKB analysis and numerical solution of the eigenvalue 
problem-that these modes are stable and have spatial scales comparable to the size of 
the disc. We show that the eigenvalue equation can be mapped to a Schrodinger-like 
equation. Analysis of this equation shows that all eigenmodes have discrete spectra. We 
find that all the models we have considered support negative frequency eigenmodes; 
however, the positive eigcnfrcquency modes are only present in power law discs, albeit 
for physically uninteresting values of the power law index j3 and barotropic index 7. 
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1 INTRODUCTION 

, Low-mass discs orbiting massive compact bodies are a fea- 

■ ture of many astronomical systems. When the dynamics of 
a disc is dominated by the Newtonian gravitational force of 

■ the central body, the disc may be considered nearly Kep- 
lerian. In a purely Keplerian potential eccentric orbits do 

, not precess because the orbital frequency is equal to the 
epicyclic frequency. In a nearly Keplerian disc there is a 
small difference between the orbital and epicyclic frequen- 
cies. This could be due to the self-gravity of the disc, ther- 
mal pressure in a gas disc, and random motions in a col- 
lisionless disc. This difference in frequencies manifests as a 
precession of eccentric orbits at rates that are small com- 
pared to the orbital and epicyclic frequencies. Then the disc 
may be able to support large-scale, slow, lopsided modes 
dKatdll983l: ISridhar. Sver fe Toumalll999l : iLee fe Goodman! 
Il999l ; ISridhar fe Toumal Il999l ). In the linear regime, these 
modes have azimuthal frequen cy, m = 1 , whos e first sys- 
tematic investigation is due to iTremainel (|200ll ) . He stud- 
ied slow modes in various types of discs (fluid, collision- 
less and softened gravity), with the focus largely on the ef- 
fect of the self-gravity of the disc. In particular, a WKB 
analysis was used to show that the fluid disc can support 
large-scale slow modes when the Mach number, ', is much 
larger than the Toomre Q parameter (both parameters are 
defined in § 2). The assumption behind this analysis is that 
the self-gravity of the disc dominates fluid pressure. How- 
ever, such is not the case for thin accretion discs around 
white dwarfs and neutron stars. Indeed, for a disc around a 



white dwarf |Frank. King fe Rainei [20021 ) . we can estimate 
^ ~ 50 and Q ~ 10 10 ; hence the analysis of (|Tremainel 
2001) is not directly applicable. The goal of this paper is to 
study large-scale m = 1 slow modes in thin accretion discs, 
where Q > M > 1. 



In § 2 we use the WKB approximation to establish 
that pressure, in the absence of self-gravity, can enable slow, 
m — 1 modes in thin accretion discs. The linear eigenvalue 
problem for slow pressure modes (or "p- modes") is formu- 
lated in § 3 for a flat barotropic disc, which is axisymmet- 
ric in its unperturbed state. When appropriate boundary 
conditions are chosen, the eigenvalue equation reduces to 
a Sturm-Liouville problem. Since the differential operator 
is self-adjoint, the eigenfrequencies are real: therefore all p- 
modes are stable, and the eigenfunctions form a complete set 
of orthogonal functions. We also map the eigenvalue equa- 
tion into a Schrodinger-like equation, which is useful in the 
interpretation of our numerical results. In § 4 we present nu- 
merical results for a variety of discs, namely (i) an approx- 
imation to the Shakura-Sunyaev thin disc, (ii) the Kuzmin 
disc which is more centrally concentrated, (iii) power-law 
discs. Of particular interest is the nature of the eigenfre- 
quency; whether it is positive or negative. This has bearing 
on the excitation of these modes, because they are stable and 
will not grow spontaneously through a non- viscous instabil- 
ity. Comparison with earlier work, summary and conclusions 
follow in 8 5. 
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2 SLOW PRESSURE MODES 

Consider a flat thin disc of fluid orbiting a central mass M. 
The fluid is assumed to be barotropic and the disc is de- 
scribed by a surface density profile E. In our analysis we 
ignore viscous forces, assuming that they adjust to main- 
tain a quasi-stationary flow with a small radial velocity, and 
have little effect on the perturbed flow. Thus, we start with 
the continuity and Euler equations in cylindrical polar co- 
ordinates 



8vr 8vr v</, 8vr v^ 2 GM d 

~~dT +VR ~dR + ~~R~!ty ~ ~R~ = W ~ dR 

dv,f, dvj, Vj, dvcf, vrv^ 

dt dR R 86 R ~ Rdd 



(* + &)> (1) 



where vr and are the radial and azimuthal components 
of the fluid velocity, h is the enthalpy per unit mass, and $ is 
the gravitational potential due to the disc. For a barotropic 
fluid with an equation of state given by p — DE 7 (where 
D > is a constant), the isentropic sound speed and en- 
thalpy are given by 



h = 



7-1 



7- 



(2) 
(3) 



2.1 Precession rate in the unperturbed disc 

We assume that the radial velocity of the unperturbed flow 
is much smaller than the azimuthal velocity and set it iden- 
tically equal to zero; this assumption is justified below at the 
end of § 2.2. The unperturbed disc is assumed to be axisym- 
metric, therefore all <f> derivatives are set to zero. Gas flows 
along circular orbits, with centrifugal balance maintained 
largely by the gravitational attraction of the central mass 
(with small but non trivial contributions from gas pressure 
and disc self-gravity) . The azimuthal and radial frequencies, 
f2 > and k > respectively, associated with nearly circular 
orbits are given by 



Q 2 



GM 1 d _ 



R dR 

d 2 



(4) 



2 GM df _ , . 3 d ., 

K = lP + dJp^ + ho)+ RdR ( ' i ' + ho) > 
where the subscript '0' indicates unperturbed quantities. 
The Mach number of the flow J((R) = Rfl(R)/c s (R) > 1. 
The dominant contribution to both Q, 2 and k 2 is due to the 
central mass, with small corrections coming from the disc 
self-gravity ($0) and enthalpy (ho)- Let us define the small 
parameter e < 1 as the larger of (E -R 2 /M) and (h Q R/GM). 
The apsides of the nearly circular orbit of a fluid element 
precesses at a rate given by, 



O - 



1 / d 2 2 d > _ 
2U\dR 2 + RdR )(*> + 



(•>) 



Note that, in contrast to lTremainel (|200lh . we have retained 
the contribution from gas pressure (i.e. enthalpy). In fact, in 
thin accretion discs around compact stars, disc self-gravity 
is negligible and the contribution to zu is almost entirely 



from gas pressure. The goal of this section is to establish 
that these discs have large-scale slow modes driven only by 
pressure. In the WKB analysis of linear modes giv en below 
we follow the presentation due to iTremainel (120011 ) . 



2.2 The WKB approximation 

We consider linear perturbations (of the velocity, surface 
density etc) of the form 



A(R) exp 



j k(R)dR + m. 



cf> — U)t 



(6) 



where k(R) and m are the radial and azimuthal wavenum- 
bers, respectively, and uj is the angular frequency of the 
mode. In the tight-winding limit where \k(R)R\ 3> 1, a 
dispersion relation between u> and k(R ) can be derived 
|Safronovlll960l ; iBinnev fc Tremaineilioog ): 



mQ) 2 — k 2 



27rGE |fe| +c 2 k 2 . 



(7) 



The disc is stable to axisymmetric (m = 0) perturbations if 
and only if 



C a K 

ttGEc 



> 1. 



(8) 



This is rea dily satisfi e d in t hin accretion discs around com- 
pact stars. ITremainel l)200ll ) showed that the dispersion re- 
lation for m = 1 modes can be written as 



uj = vj + 



TrGEojfcj 



In ITremainel (|200lf ) it is argued that, when the pres- 
sure is negligible compared to disc self-gravity (i.e. c s ~ 0), 
the WKB dispersion relation of equation © admits large- 
scale '\k(R)R\ ~ 1' modes with frequencies u> ~ w ~ 
(Eo-R 2 /M)fi. It may be verified that the condition of negli- 
gible pressure implies that the Mach number ^£ 3> Q- How- 
ever, as we have argued in the introduction, this inequal- 
ity is violated for thin accretion discs around compact stars 
where the opposite is true, i.e. Q 3> 3> 1. Hence we need 
to co nsider a situation that is complementary to the analy- 
sis of ITremainel (|200lh . Disc self-gravity being negligible in 
accretion discs, the precession rate is determined entirely by 
the gas pressure. Then equation <(Sj can be written as 



1 



d 2 h 2 dh 
dR 2 + R~dR 



0(1/^ 4 



n 



(10) 



fsee lKatol l|l983h ) and we can approximate equation ((9]) as 
c 2 k 2 1 



-0(w 2 ,u} 2 ) 



(11) 



c s k 

~2n 1 n" 

For a disc with a non zero inner radius, eigenmodes must 
satisfy the Bohr-Sommerfeld quantization condition, given 
by 



jk{R)dR= (n+ ^ 2tt, 



n = 0,1,2,. 



(12) 



Thus, there exists a prima facie case for modes with fre- 
quencies uj ~ uj, with radial wavenumbers, 

1/2 

k{R) 



1 

R 



(13) 



comparable to the radial scale of the disc. However, this 
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tentative conclusion is based on a WKB analysis which may 
not be valid for modes with k(R)R ~ 1. This is the mo- 
tivation for our studies of the eigenvalue problem for slow 
modes given below. Henceforth we ignore disc self-gravity 
altogether and consider only the effect of gas pressure. 

We now justify the assumption made at the begin- 
ning of § 2.1, that the radial velocity in the unperturbed 
disc is small, and may be ignored when studying slow 
modes. The time scale of radial spreading of the disc is, 
tvis ~ R 2 /u ~ ^ 2 /(aQ), where a <C 1 is the Shakura- 
Sunyaev viscosity parameter. The frequency of a slow mode 
is W ~ w ~ fL#~ 2 . Therefore ut V i s ~ a 1 3> 1 implies that 
the radial spreading occurs over many slow mode periods. 



3 FORMULATION OF THE EIGENVALUE 
PROBLEM 

3.1 Eigenvalue equation 

The linearized Euler, continuity and enthalpy equations that 
govern the perturbed flow are 



dvr 



dt 
dt 



- 2Q(R)v^ - 
2B(R)v R1 



~dT "R dd> 



dhi 
' OR ' 
1 dhi 
~R dcf> ' 

l.JL moVRl) 



In 



2 El 
: c s0'^ _ ; 



(14) 
(15) 
(16) 
(17) 



where the subscript '0' stands for the unperturbed quantities 
and T' for the first order perturbed quantities. The Oort's 
parameter B(R) is related to the epicyclic frequency through 
k 2 (R) = -4S}(R)B(R), and d/dt = (d/ffi + fld/d(t>) is the 
convective derivative with respect to the unperturbed flow. 

We consider non-axisymmetric perturbations with az- 
imuthal wave number m = 1, of the form T\ = 
T a (K) exp [i(<f> — Lot)], where 2i stands for any perturbed 
quantity. Substituting this form in equations (f!4|) . 1(15) , (fl6)) 
and |(T7) yields 



l(Q — Lu)v Ra — 2flV4, a + = , 

dR 



2Bv Ra + % = 
R 



., . 2E0 1 d . , 

- Uj)T, a + —j^-V^a + — — (RT.oVRa) 

h — 2 2 ~-' a 



(18) 
(19) 
(20) 
(21) 



Solving equations (|18[1 and p9|l for the velocity amplitudes 
we obtain 



VRa 



I 

A 



1 

A 



. .dh a 20. 
^-^iR+ll^ 

dh a Fl-u 
- 2B m + -FT 11 " 



(22) 
(23) 



where A 



- (fi — uj) 2 . These equations, along with equa- 



tion (|21[) . when substituted in (|20p yields 



d 2 
dR? 



+ 



d i 



i?E 



_d_ 

dR 



29. 



— h 

R(Q-u) \ dR n 



E fi 



+ 
1 

R 1 



ha = 



h a A 



(24) 



which is the eigenvalue problem for undriven modes, with 
eigenvalue uj and eigenfuncti on h a . This equation i s a spe - 
cial case of equatio n (13) of iGoldreich fe Tremaine] (1 19791 ); 
iTsang fc Lail (2009). where m = 1 , and the external and self 
gravity perturbations are set equal to zero. It can be noted 
that this equation becomes singular at SI = lo and A = 0. 
The former corresponds to the corotation resonance and the 
latter corresponds to the Lindblad Resonances (LR) . Below 
we discuss the validity of equation (|24|l at LR; a similar anal- 
ysis holds for the singularity at corotation radius but we do 
not discuss it in this paper since, as is argued later, for the 
slow modes the corotation radius has to lie outside the disc. 

The system of equations (|18jl — (|2ip describe an un- 
driven, autonomous system. At the Lindblad resonance the 
algebraic equations (|18|) and (|19p . become indeterminate if 
no conditions are imposed on the enthalpy perturbations ha- 
lt is easily seen that these equations become consistent if 



i(Q -u) _ 2iQ, _ dhajdR 
2B ~ (Q-u) ~ ih a /R 



(25) 



The first equality follows from A = 0. Rearranging the sec- 
ond equality yields 



_d_ 

dR 



(R 2 h a ) 



R 



dh a 
~dR 



0. 



(26) 



This condition must be satisfied at the Lindblad resonances 
for all undriven modes. However, equation (|26[) may not be 
satisfied if the disc is driven by external forcing and may lead 
to curious dynamics around the LR and transport of angu- 
lar momentum away from th e LR due to external torquing 
IIGoldreich fc Tremaine! Il979h . In this work we confine our 
investigations to free modes of an undriven disc, therefore, 
as the above discussion shows, nothing special happens at 
the LR. 

Slow Mode Approximation: We now make the 
ansatz that the perturbed flow supports frequencies that are 
small in comparison to the circular frequency, i.e., \uj\ <g; Q,. 
Therefore, when u / 0, the disc must be finite, with outer 
radius such that the orbital frequency at the outer edge is 
much greater than |w|. Applying the slow mode approxima- 
tion to equation (|24[) we obtain 



2 30 R 3/2 d 



E n 



E « dR \ R 3 / 2 A dR 



d (R 2 h a 



= R^ha 



(27) 



Similar to equation (124|l , equation (|27[) too is singular at the 
Lindblad resonances, however the singularity at the coro- 
tation radius has gone away since this equation has been 
derived under the slow mode condition, \uj\ <<£ Q,. The con- 
dition A = implies that at some radius, either ui — fl — k 
or lo — Q + k,. Since n ~ Q, + 0(e 2 ), we see that the second 
equality cannot be satisfied under the slow mode approxi- 
mation. It is straightforward to see that the radius where 
this would be satisfied would be larger than the corotation 
radius due to the fact that the Keplerian circular frequency 
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falls off monotonically with radius. Therefore, there are no 
outer Lindblad resonance singularities for slow modes. How- 
ever, the Inner Lindblad Radius (ILR), where u> = tb(R), 
could very well lie inside the disc. Due to the fact that the 
disc surface density is completely arbitrary, there could in 
general be more than one ILRs. To make the problem well 
posed under the slow mode approximation, at the ILRs, the 
condition (1261) reduces to 



= o. 



(28) 



We shall see later that the numerical solutions of equation 
(|30[) satisfy this condition and the velocity amplitude at the 
ILR remains finite; thus the linear approximation remains 
valid and nothing special happens at the ILR. 

From equation (|10[) . zu/fl ~ 0(^~ 2 ). This allows us to 
approximate k 2 ~ Q 2 , leading to B ~ —$1/4. Using these in 
equation (|23|) we obtain 



v <t>a — r, A 



Q_ 
2~A 



dh a 
~d~R 
1 d 



2hg 

R 



(R 2 h a ) 



2A R? dR v " ' (29) 

Differentiating equation (|27|l and using equation (|29[) we ob- 
tain, 



_d_ 

dR 



2 j-,3/2 

So SI 



d& 
dR 



2R?I 2 , 



■ w)e = o, 



where we have used the variable O = R 1 ^ 2 T,ov ( t >a and 



(30) 



A ~ 2fi(w - 
proximation. 



zo), which is valid under the slow mode ap- 




Radial coordinate (x) 



Figure 1. The effective potentials for our barotropic approxima- 
tions to the SS disc. The solid line corresponds to Vi(x) and the 
dashed line to V2(ai), described in § 4. The positive values of A 
can provide both Type I and Type II eigenvalues as described in 
the text. The negative values, although seemingly allowing eigen- 
states lead to no such solution. Note that positive values of A 
correspond to the negative frequency modes. 



3.2 Slow modes as a Sturm-Liouville problem 

Before we proceed to specific examples, we have to choose 
the boundary conditions that we impose to solve equa- 
tion (|30[) . We first cast the equation in a dimensionless 
form by choosing a radius R+, at which we evaluate var- 
ious quantities, , c s * , to* and $7*. We introduce the pa- 
rameter x — R/Ri,, and similarly for a quantity H, we use 
H' = H/ H^, leading to the Sturm-Liouville form of the eigen 
equation 



^(P(x)^) Kilr) -AWl ,■))(-) = (I. 



(31) 



where 



P(x) = 



Q(x) = 



c' V /2 



W(x) = 



2x 3 ' 2 



2x 



3/2 T 



e 



1/2 



v<f, a T,' . 



In equation (|31|) . A = — u;„#* 2 /f2* is defined with a neg- 
ative sign to make an explicit correspondence with the 
Schrodinger's equation, to be introduced in § 3.3. Hence- 
forth we reserve the term "eigenvalue" for A, and use either 
"frequency" or "eigenfrequency" for u. 

We consider discs with an inner edge at 7?i n ncr, and an 
outer edge at i? ou tor, and we choose i?* = iiinncr- We now 
argue that the parameters of the disc and the central mass 
for astrophysically interesting discs are such that the slow 
mode condition can be easily satisfied everywhere inside the 



disc. In a Keplerian disc the slow mode condition \u)\ <sK $1 
is satisfied everywhere in the disc if it is satisfied at the disc 
outer radius, this leads to 



Ro 



R, 



3/2 



|A| 



(32) 



where we have used A = — u^ + 2 /n*, and fi(-Router) = 
f2*(i?*/-R ou tcr) 3 ^ 2 . Typical expected values for are in 
the range 10 4 -10 6 (Frank. King fc Raind |2002| ) . Most ex- 
amples we consider have surface densities that decline by 
Router /R* — 30-50, therefore we see that for an eigenmode 
to be slow through out the disc, |A| has to be much smaller 
than ~ 10 5 -10 9 . We shall see in the examples that this con- 
dition is comfortably satisfied. 

We integrate the eigen equation in the range 1 < x < 
Pouter- We assume that the perturbations obey the boundary 
conditions, 



9(1) = 9(a;oute r ) = 0. 



(33) 



We note that these boundary conditions make the differ- 
ential operator in equation (|31[) self-adjoint: therefore the 
eigenvalues A are real, and all slow p-modes are stable. The 
complete set of eigenfunctions also form a complete basis; 
however, we note that not all eigenvalues are slow, and thus 
we do not expect this set to describe the evolution of arbi- 
trary perturbations, but only the ones that obey the slow 
mode condition, $1 ^> |w|. 
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3.3 Effective potential and WKB approximation 

In the usual WKB approximation we substitute the trial 
solution 



Q(x) = A(x) exp 
in the following equation: 



— / kdx 



(34) 



fi 2 -^{P(x)^j+(Q(x) + \W(x))e = 0. (35) 

Here p is an ordering parameter which is finally set equal to 
unity. A(x) and k(x) are the amplitude and the wavevector 
respectively. Collecting terms of zeroth order in [i leads to 
the dispersion relation 



k 



Q(x) + XW(x) 
P(x) 



2Q 



(vj - 



(36) 



which is identical to equation However, we find that 

this dispersion relation, together with the Bohr-Sommerfeld 
quantization condition of equation (I12|l predicts eigenval- 
ues that compare poorly with those obtained from nu- 
merical integration of the Sturm-Liouville equation. Hence 
we have reformulated equation (|3ip . using new variables 
rj{x) = ^/P(x) and * = y / P(x)Q. Then equation J3TJ takes 
the Schrodinger-like form 



where 



K 2 {x) 



r] 2 (x) 



*" + K 2 (x)^ = 0, (37) 



\Q(x) + XW(x)-r l (x)v' (x)] , (38) 



which on defining V(x) = (— Q(x) 
be written as 



r,(x)r,"(x))/W(x) 



K 2 (x) 



W(x) 
r/ 2 (x) 



[A-V(z)] . 



(39) 



This dispersion relation differs from equation |TT} and seems 
to better describe the numerical solutions, giving a match 
with the numerically obtained eigenvalues up to a few per 
cent, as can be seen in Table 1. 

Note that K 2 (x) in equation (|39|l differs from the 
standard form for the Schrodinger equation by the factor 
W(x)/rj 2 (x). However, it is very useful for discussions of 
the turning points, where K 2 (x) = 0, separating classically 
accessible regions from the forbidden ones. The solution is 
oscillatory where K 2 > 0, implying A > V(x), and non- 
oscillatory otherwise. Since the disc is finite the eigen spec- 
trum is always discrete and there are two distinct types of 
spectra: 

Type I: This occurs when there is at least one turning point 
within the disc. In the case of a single turning point we 
can have oscillatory solution on either side of the turning 
point, depending on the form of K 2 (x). If there are more 
than one turning point then we could either have oscillatory 
behaviour confined between the turning points or outside, 
such as the case of the SS disc discussed below. 
Type II: This occurs when there are no turning points 
within the disc and the discreteness of the spectra depends 
entirely on the size of the disc. 

To obtain real eigenvalues we need to consider the possi- 
bility of satisfying K 2 > in a bounded region, which could 
either be bounded by one or both of the disc boundaries 



d 
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20 
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Figure 2. The effective potential for the Kuzmin disc. Positive 
values of A lead to a discrete spectrum of both Type I and Type II 
eigenvalues. There are no negative eigenvalues. 



or by turning points. Thus a useful first step is to plot this 
potential for the problem at hand. If the potential allows 
regions that can support bound states, we search for solu- 
tions numerically, and to verify our results we use the WKB 
approximation. 

Let us consider the case of the SS disc shown in Fig [1] 
On the negative side the potential blows up at the inner edge 
at x — 1, where the perturbations are assumed to vanish. If 
A < is to be a valid eigenvalue, then it must satisfy the 
quantization condition 



i; 



K(x) dx — 



n = 0, 1,2, 



(40) 



where, x — a > 1 is a turning point which separates a classi- 
cally accessible region (I) on the left from a forbidden region 
(II) on the right. The case of positive eigenvalues, A > 0, is 
more interesting. If the eigenvalue is positive and smaller 
than the maximum value of V(:r), then there are two turn- 
ing points (say, x — a,b with b > a) separating three dis- 
tinct regions. The classically forbidden region lies between 
x — a and x = b, separating the two classically accessible 
regions (l,o) and (6, Xouter)- If A is greater than the maxi- 
mum of V(x), then whole disc is classically allowed. Below 
we present numerical results on eigenvalues and eigenfunc- 
tions, and use WKB approximation to understand them. It 
turns out that WKB approximation is very useful even for 
the case of small quantum numbers. 



4 NUMERICAL RESULTS 

As discussed in the last section, we consider finite disc and 
we expect the spectra of equation (|31|l to be discrete. These 
modes would have observational consequences since they 
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Figure 3. The allowed region (hatched) in the /3-7 space to ob- 
tain Type I negative eigenvalues (corresponding to the positive 
frequency modes) for the power law discs (case (ii)). Most of the 
region is either unphysical or uninteresting. 



would rotate at a definite frequency around the disc. The 
perturbations in the enthalpy would lead to azimuthal vari- 
ations in the temperature and density across the disc, which 
might be observable depending on the amplitude of pertur- 
bations. 

Since very little is known about the surface density pro- 
files of the discs we carry out a simplistic calculation based 
on certain standard forms of the disc as test cases. We also 
consider the generic power law profile. Some of the profiles 
considered below are formally infinite in size, however as 
noted above, we need to keep the disc finite. This would im- 
ply that the surface density would abruptly fall to zero at 
the outer disc radius. This is unphysical and we expect that 
there would be a thin transition region that would deviate 
from the density profile being considered near a; ou tcr. If the 
eigenvalues and eigenmodes are not very sensitive to x on t eT , 
then this is not an issue and in our numerical investigation 
we indeed find this to be true. In all the examples below 
Pouter = 50. In Table 1 we give the first a few eigenvalues for 
Xoutei = 50 and Xoutei = 200, as A50 and A200 respectively. It 
can be noticed that eigenvalues do not change substantially. 

Shakura-Sunyaev (SS) discs: We first consider 
the standard model of an accretion disc pr oposed by 
Shakura and Sunyaev (|Shakura fc Suny acv 1973). The sur- 
face density and temperatu re of this disc are given by (see 
iFrank. King fc Raind (|2002h l: 



(41) 



/ = 



Rio 



1/2 



1/4 



(43) 



Mi6 is the mass accretion rate in the units of 10 g s , mi 
is the mass of disc in solar mass units, Rio is radius in the 
units of 10 10 cm, R* is the radius of the central object in 
the units of 10 10 cm, and a is the Shakura-Sunyaev viscosity 
parameter. 

Although the SS disc is not based on a barotropic 
model, we find that a barotropic disc with index 7 = 2 
serves as a reasonable approximation. We have considered 
two cases: 

(i) Choosing Eo = Ess of equation (I4ip . and deriving 
the temperature profile for 7 = 2, we find that the effective 
potential is 

344 - 590^/x + 225x 



Vi(x) 



800 :r 9 / 4 (l 



(44) 



(ii) Choosing To = Tss of equation (|42l) . and deriving the 
surface density profile for 7 = 2, we find that the effective 
potential is 



Va(z) 



12 -22 ^ + 9a; 
32 a; 9 / 4 (l- 



(45) 



where we have used the natural length scale given by the 
inner disc radius, which we have adopted for conversion of 
our eigen equation into a dimensionless from. As Fig {1} 
shows, Vi(x) and Va(a;) are quite similar to each other. The 
effective potential blows up at the inner edge of the disc and 
steadily climbs up above zero and then decreases asymptot- 
ically. In the immediate vicinity of the inner edge a negative 
A gives an oscillatory solution, and to the right of it is a 
classical turning point. Beyond the turning point the solu- 
tion is exponentially decaying. This suggests that discrete, 
Type I, negative eigenvalues might exist, however, both the 
numerical search and the WKB quantitation condition (|40fl . 
fail to find such discrete eigenvalues. For < A < 0.01, we 
find discrete, Type I eigenvalues for which the oscillatory be- 
haviour is outside the region bounded by two turning points. 
For A > 0.01, we find discrete eigenvalues of Type II, where 
the separation between neighbouring eigenvalues decreases 
with increasing values of the outer radius of the disc. 

The Kuzmin Disc: In contrast to the SS discs, the 
Kuzmin disc has a centrally concentrated surface density, 
and hence offers a distinct case in which to study slow 
modes. The surface density profile in this case is given by 

S <*>= M^+V ' (46) 

where Md is the mass in the disc and a is the core radius. 
The surface density extends all the way to R = 0. If we 
take a as the size of the inner radius of the disc then we 
can rescale our equation by choosing R+ = a, leading to the 
effective potential 



V(x) = 



3(z 4 (5 - 4 7 + 3 7 2 ) + x 2 {6 - 87) + 1) 



(47) 



/ K , (42) 



8x 1 / 2 {l+x 2 )( 3 ->+ 1 ^ 2 

In Fig (fj]) we have plotted this potential for 7 = 4/3. As 
may be seen, this case qualitatively resembles Fig {l}, and 
the discussion in § 2.4 carries through. For negative values of 
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Figure 4. The allowed region (cross-hatched) in the /3-7 space 
to obtain discrete positive eigenvalues (corresponding to the neg- 
ative frequency modes), corresponding to Type I eigenvalues for 
case (iii), for the power law discs. 

A, we can infer from Fig ((2} that K 2 (x) < 0, so wave like so- 
lutions are not possible. For positive values of A in the range 
< A < 0.3, we can have Type I eigenvalues, but this region 
of A is further divided into two parts: for < A ^ V(l), 
there is only one turning point and oscillatory behaviour is 
possible for radius greater than the turning point, and for 
V(l) < A < 0.3, there are two turning points and oscillatory 
behaviour is possible outside the region bounded by the two 
turning points. For A > 0.3 the eigenvalues are of Type II. 
This behaviour is confirmed by numerical integration of the 
eigenvalue equation. We can also admit values of 7 other 
than 4/3. However, the effective potential in equation (|47p 
retains the general shape of Fig and the conclusions 
stated above remain valid. 

Power Law D iscs: Certain physical models (e.g. 
iNaravan fc Yil l|l994l )) naturally lead to scale invariant discs 
that follow a power law profile. Although in these models 
there is no associated length scale, we choose to truncate 
the disc at the inner disc radius thus leading to the form 



E = £* x p , 



(48) 



where, E* = E(i?*) and /3 is the power law index. The po- 
tential is given by 



where 

C(/?,7) 



V(x) = C{p,i)x v , 



(3 + (-4 + 7(4 + 70))) 



(49) 

(50) 
(51) 



7/3 - - 1/2 . 
There are four distinct possibilities: 
i) C(/3, 7) < 0, v(fi, 7) < 0: The region in /3-*f space 
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Figure 5. Two eigenfunctions (Type-I eigenvalue) obtained for 
the power law case (case iii) with f) = 7 = 2. The higher quantum 
number leads to a more radially extended eigenmode. Star marks 
the position of the turning point. 
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Figure 6. The allowed region (hatched) in the /3-7 space for the 
positive Type I eigenvalues (corresponding to negative frequency 
eigenmodes), for case (iv) for the power law discs. 
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n 


WKB A 


Numerical A50 


Numerical A200 





13.42 


13.86 


13.86 


1 


30.15 


30.73 


30.73 


2 


52.31 


53.06 


53.06 


3 


79.99 


80.92 


80.91 


4 


113.22 


114.32 


114.31 


5 


151.99 


153.29 


153.26 


6 


196.31 


197.83 


197.75 


7 


246.18 


247.96 


247.80 


8 


301.60 


303.72 


303.40 


9 


362.57 


365.16 


364.54 


10 


429.09 


432.31 


431.24 



Table 1. The eigenvalues for the power law model with /3 = 7 = 2 
are tabulated here for comparison between those obtained nu- 
merically and those obtained using the WKB approximation. 
The columns A50 and A200 are the eigenvalues corresponding to 
Pouter = 50 and 200 respectively. The match between numerical 
and WKB eigenvalues is within a few per cent, and remains good 
even for small quantum numbers. The eigenvalues are not very 
sensitive to precise value of £ utcr 



tended, with well separated peaks. For A > V(l) there is no 
turning point and the eigenvalues are of Type II. 



satisfying these conditions is plotted in Fig ((3j). The regions 
where these conditions are satisfied are those with (1) f3 
positive and 7 negative or (2) very small negative j3 and 
very large 7. Both these cases are unphysical. Note that 
there is a region with large, positive /3, and small, positive 
7, and these are also physically uninteresting. 

(ii) C(J3, 7) < 0, u(J3, 7) > 0: It can be verified that for 
v{j3, 7) > C(/3, 7) is always positive and hence it is impos- 
sible to satisfy these conditions. 

(iii) C(/3, 7) > 0, v(/3, 7) > 0: These two constraints give 
us a region in the space, displayed in Fig Q , that admits 
physically reasonable values. For V(:ri nncr ) < A < V(x outC r), 
we have one turning point admitting Type I eigenvalues. In 
Fig ((5]) we plot two examples of eigenmodes for the case 
= 7 = 2. The eigenfunctions for small quantum numbers 
are found to be centrally concentrated while they extend to 
larger radii for larger quantum numbers. It should be noted 
that our solutions are regular at the turning pointfl The 
first few eigenvalues for this case are tabulated in Table 1, 
where we find an excellent match with the WKB eigenvalues. 
Outside this range of positive A there are no turning points 
and hence only Type II eigenvalues are possible. 

(iv) C(/3, 7) > 0, 7) < 0: The region in /3-j space sat- 
isfying these constraints is plotted in Fig |6j. For < A < 
V(l), we have one turning point admitting Type I eigenval- 
ues. Here the classically accessible region is bounded by the 
turning point on left and x ou t er on right. Numerical solution 
for P — —2 and 7 = 2 are plotted in Fig @. Comparison 
with Fig (0 shows that the solutions are more radially ex- 



The positions of the ILRs, obtained from, A = — Q(x)/W(x) 
are slightly different from the turning points obtained from the 
equation A = -(Q(x) - T)(x)t)" (x))/W(x) . 



5 DISCUSSION AND CONCLUSIONS 

We have presented a theory of slow m = 1 linear pres- 
sure modes (or "p-modes") in thin accretion discs around 
massive compact objects, such as white dwarfs and neutron 
star. These modes are enabled by the small deviation from 
a purely Keplerian flow, due to fluid pressure rather than 
disc self-gravity. For simplicity we have taken the fluid to be 
barotr opic. Our formulation largely follows that of lTremainel 
althoug h there is a key difference: using the WKB 
approximation, ITremainel l|200lf) argued that fluid discs for 
which disc self-gravity dominates fluid pressure can support 
slow modes, if the Mach number is much larger than the 
Toomre Q parameter. This condition may be satisfied in 
relatively cool discs, but not for thin accretion discs around 
white dwarfs or ne utron star s . In t hese discs, Q S> ^ ^> 1, 
and the analysis in ITremainel (2001) does not apply, because 
disc self-gravity is negligible when compared with fluid pres- 
sure in thin accretion discs. This implies that the precession 
rate of the apsides {zu) of a fluid element in a nearly cir- 
cular orbit is determined by the fluid pressure; to order of 
magnitude, vj ~ {Q,/^ 2 ). 

A WKB analysis was used first to argue that thin ac- 
cretion discs can support large-scale (k(R)R ~ 1), m = 1 
p-modes wit h small angular frequencies, lo ~ vj ~ (fi/„# 2 ). 
As noted by ITremainel l|200ll ) . these long wavelength modes 
may dominate the appearance of the disc, and are not ex- 
pected to be damped by viscosity. We derived an eigen equa- 
tion for slow linear modes and showed that it is identical to 
a Sturm-Liouville problem. The differential operator being 
self-adjoint implies that the eigenvalues are all real, so that 
all slow p-modes are stable. This corresponds to the result in 
ITremainel (|200ll ) that all slow modes of the softened gravity 
disc are stable. We solved the Sturm-Liouville problem nu- 
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merically for a variety of unperturbed discs, and summarise 
our results below. 

(i) The first corresponds to two different kinds of 
barotropic approximations to the Shakura-Sunyaev thin 
disc, which have modes with negative eigenfrequencies. 

(ii) The second is the Kuzmin disc, which is more cen- 
trally concentrated. As earlier, this too supports only nega- 
tive eigenfrequencies. 

(iii) Power-law discs can support modes with negative 
eigenfrequencies for reasonable values of (5 and 7. For cer- 
tain combinations of these parameters power law discs can 
support positive eigenfrequencies as well; however the range 
of parameters turn out to be physically uninteresting. 

If slow modes are stable, it is necessary to consider how 
they could be excited. Since they have azimuthal wavenum- 
ber m — 1, we need look for excitation mechanisms which 
possess the same symmetry, at least in the linear limit. One 
possibility is from the stream of matter from the secondary 
star that feeds the accretion disc. When viewed in the rest 
frame of the primary, the region where the stream meets 
the outer edge of the disc rotates in a prograde sense with 
angular frequency equal to the orbital frequency of the bi- 
nary system. Since the orbital frequency of the binary can be 
much smaller than the orbital frequency of the gas in the ac- 
cretion disc, there is the possibility of the resonant excitation 
of a slow mode if it has positive frequency. In the cases we 
have considered, we find only negative frequencies belong- 
ing to a discrete spectrum, allowing only for non-resonant 
driving. 

There is, however, another alt ernative that does not rely 
on external sources of excitation. IZhang fc Lovelace] (|2005l ) 
have studied linear waves in thin accretion discs and ap- 
plied their theory to slow m — 1 modes around black holes. 
They used the pseu d o-New tonian gravitational potential of 
IPaczvnskv fc Wiital l|l98d ) to model the general relativis- 
tic effects due to a Schwarzschild black hole. In this case 
■uj is due to the deviation of the pseudo-Newtonian poten- 
tial from a Kepler potential, so their slow modes are not 
driven by pressure (as is true in all the cases we have con- 
sidered); hence more detailed comparison with our work is 
not possible. What is interesting is that they find that their 
slow modes to have negative energy and angular momen- 
tum, and suggest that slow modes may be excited sponta- 
neously through the action of viscous forces. This possibil- 
ity should be examined in the context of the p-modes we 
have studied. However, this would require reformulating the 
eigenvalue problem taking into account viscous effects. 
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